Long-living Bloch oscillations of matter waves in optical lattices 
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It is shown that by properly designing the spatial dependence of the nonlinearity it is possible to 
induce long-living Bloch oscillations of a localized wavepacket in a periodic potential. The results 
are supported both by analytical and numerical investigations and are interpreted in terms of matter 
wave dynamics displaying dozens of oscillation periods without any visible distortion of the wave 
packet. 
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The phenomenon of Bloch oscillations (BO), predicted 
by Bloch in 1928 in his celebrated paper on the dynamics 
of a band electron in a steady electrical field repre- 
sents a problem of non exhausted interest. Besides solid 
state physics, where the phenomenon has been observed 
only in recent times after the development of the super- 
lattice technology it is now possible to observe BO 
also in other fields such as nonlinear optics, using light 
beams in arrays of waveguides [sl or in photorefractive 
crystals and atomic physics, using Bose-Einstein con- 
densates (BEG) loaded in optical lattices (OLs) [1, 0, [3|- 
Apart its fundamental significance, the interest in BO 
arises mainly from perspectives of their practical appli- 
cations. In this context we mention the use of BO for 
metrological tasks, including relatively precise definition 
of h/m [8| and measurement of forces at the micrometer 
scale like the Casimir-Polder force [§) and the gravity ■ 
Recently, BO were also suggested as a tool for control- 
ling light in coupled-resonator optical waveguides llj . In 
these physical contexts BO have been observed mainly 
in the linear (as first proposed by Bloch) or quasi linear 
regimes (where the nonlinearity introduces quantitative 
but not yet qualitative changes). The present experi- 
mental settings (both in nonlinear optics and in BECs), 
however, allow for attaining essentially nonlinear regimes 
quite easily, this actually being a desirable condition for 
applications involving localized states. 

The fact that BO can exist in presence of interac- 
tions (nonlinearity) was first recognized in the context of 
nonlinear discrete systems 12[. For periodic continuous 



models of the nonlinear Schrodinger (NLS) type, such 
as ones describing matter waves in OLs, the existence 
of BO becomes more problematic because of nonlinear- 
ity induced instabilities in the underlying linear system. 
These instabilities can be simply understood by observ- 
ing that in a usual OL at the edges of an allowed band 
the effective mass has always opposite signs. This im- 
plies that if a Bloch state is modulationally stable (in 



presence of a constant nonlinearity) at one edge of the 
band, it must be necessarily unstable at the other band 
edge [l3]- Nonlinearity induced instabilities have been 
extensively investigated both theoretically [l^, [l3| and 
experimentally [15]) and have been identified as primary 
cause for the short life-times of BO of matter waves (the 
oscillation survives only a few cycles) observed both in 
numerical [T^ and in real experiments [ist . 

In this Letter we show that by properly controlling 
the instabilities of the system it is possible to achieve 
long-living BO of matter waves in OLs. In our approach 
BO are implemented using solitonic propagation (the 
situation similar to one explored in truly discrete sys- 
tems [12|). To this end we use spatial periodic modula- 
tions of the scattering length via optically induced Fesh- 
bach resonances, in order to change the stability proper- 
ties of the Bloch states at the edges of the band. These 
modulations correspond to a nonlinear lattice whose am- 
plitude, considered as a free parameter, can be used to 
eliminate instabilities from the band (similar stabilizing 
properties of nonlinear lattice were used to achieve de- 
localizing transitions in one dimensional (ID) case [l3])- 
We show that the regions of parameter space for which 
Bloch states become unstable (almost) in the whole band 
coincide with those for which long-lived BO of matter 
waves become possible, this showing the validity of our 
approach. 

We start by considering the following ID Gross- 
Pitaevskii equation 



-ipxx + ixxp + u{x)tp + gix)\iijfij, (1) 



describing an array of BEGs in linear and nonlinear tt- 
periodic optical lattices of the form U{x) — U{x -\- n), 
Qi^) = CJ(a; + tt), respectively, with the external linear 
force 7 arising from an uniform acceleration of the lat- 
tices (a similar approach is valid also for the multidi- 
mensional case). The consideration will be restricted to 
symmetric potentials U{x) — lA{—x). Since BO exist in 
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the linear regime it is appropriate to start our consid- 
erations from the underlining linear eigenvalue problem 
—(pipq/dx^ + U{x)(pq = e{q)(pq, where tpq is the stan- 
dard Bloch function corresponding to the wave vector 
q. The periodicity of the linear OL introduces a band 
structure in the spectrum e.g. the existence of chemical 
potential functions £^{9) which are periodic in reciprocal 
space q G [—1,1], with n denoting the band index. In 
the following we shall consider only the lowest chemical 
potential band, a situation typical for most BEC appli- 
cations, and omit the band index. In presence of nonlin- 
earity, stationary localized states (gap-solitons) can exist 
only if chemical potentials are inside gaps (see e.g. [l8|). 
This is not the case of a nonstationary (moving) soli- 
ton whose velocity, in the leading order, coincides with 
the group velocity of the carrier Bloch wave and mathe- 
matically can be described by a multiple-scale expansion 
provided the linear force is weak enough 7^/^ <C 1. To 
this regard we search for localized solutions of Eq. ([T]) in 
the form ip = +7V'2 + • • • , where at each order ipj 

denotes a function of a set of scaled temporal t„ = 7"/^^ 
and spatial x„ = 7"/^a: (n = 1, 2, ...), variables. Since for 
BO the solution, thought as a wavepacket of Bloch states 
ipq with the carrier wave-vector depending on the slow 
time q = q{^t), must scan the whole band, the ansatz for 
the first order ipi is chosen, in analogy with the Houston 
functions [3] of the underlying linear theory, of the form 

^i=A(T,Oe'^^*V,M(x) (2) 

where the phase £{t) is determined by the equation 
d£/dt — e{q(jt)) and we have introduced specific nota- 
tions for the slow variables t = jt and ^ = ^y^{x — v{T)t). 
This peculiarity of the asymptotic expansion will man- 
ifest at the third order of the small parameter, i.e. 
0(7'^/^), lower orders of the expansion being obtained 
in standard manner (see e.g. [13|). Since the temporal 
dependence of the wavevector is not specified, we impose 
the constraint that q must follow the well known semi- 
classical equation for the linear BO which in our scaled 
variables acquires the simplest form q = t[21|. Then, 
in the order 7'^/^ we arrive at the NLS equation for the 
slowly varying amplitude 

'^^ + 2F(^^««"^('^^'^''^^'^- 

Here we defined the group velocity v{q) — de{q)/dq, the 
effective mass M{q) ~ {2d'^e{q) / dq^)~^ and the effec- 
tive nonlinearity xil) = /o Q{x)Wq{T){x)\'^dx . We em- 
phasize that although we have indicated q as an argu- 
ment in the above definitions, the group velocity, the 
effective mass and the effective nonlinearity are func- 
tions of the slow time r. Introducing and "— " 
subindexes to denote properties at top and bottom lim- 
its of the band, we have that ip^~'^ = V5q=0i </5*-^^ = </?g=i, 
;,(^(±) = g{x)\ip^^\x)\^dx , and T^'^' > 0. As it is 



clear, if Q{x) =const, then x''~^^X^~^ > ^md from Eq.Q 
it follows that while envelope solitons are available at one 
edge of the band (the edge for which Mx < 0), they can- 
not exist at the other edge. Since for BO the soliton must 
travel along the whole band, it will necessarily reach the 
edge where it undergoes strong dispersion (due to defo- 
cusing action of nonlinearity {Mx > 0), this leading to 
destruction of BO J^] . 

From this analysis it is clear that an essential control 
parameter of the problem is cr{q) — —M{q)x{q) and 
in order to have long-living BO we must require that 
(y{q) > for all q, this assuring a focusing nonlinearity 
for the soliton dynamics along the whole band. The 
above condition can be achieved by means of a proper 
design of the nonlinear lattice and will be optimized if 
the condition a =const> is satisfied. Analytically it 
is easy to consider the general case in which (t{t) is a 
slowly varying function \da/dq\ <^ a. Indeed, introduc- 
ing the new temporal variable t = JJ" dr' /M (r') and the 
new function A = A/ ^/a we obtain the NLS equation 
with a dissipative term lAr + \A^i_ + \A\'^A — i§^A 
whose approximate solution can be found by means 
of soliton perturbation theory For the partic- 

ular case of a static soliton in the frame moving 
with the velocity v{t), the solution reads A — 

Ao cxp (r(?) + ^Al /; ^M^dr) sech (AoHr)]^/^^) 
with ^0 the constant amplitude of the soliton. 

In order to check the above predictions we have per- 
formed direct numerical simulations of Eq. ([T]) with 
U{x) — — Fcos(2x) and g{x) = Gq + GiC0s(2a;), where 
V and Gi denote the amplitudes of the linear and non- 
linear lattices, respectively, Gq is the "average" nonlin- 
earity, which in all numerical simulations reported be- 
low is chosen to be Go = —0.777. We denote by G^^^ 
the value of Gi at which effective nonlinearity x''^^ be- 
comes zero. Respectively if Gi ^ G^^'' we have that 
X^^-* ^ and the condition for the existence of enve- 
lope solitons (alternatively, the instability of the Bloch 
waves) at the both band edges M'^^^x^^'' < is met 
when g\^^ < Gi < g[ \ This corresponds to the do- 
main of parameters between the lines in FigUwhere G^^' 
vs V are depicted (e.g. point B in FiglT]). Consequently, 
domains of parameters Gi < g\^'' (e.g. point A in FiglT|) 
and Gi > G^""* (e.g. point C in Fig[T|) allow for the exis- 
tence of small amplitude envelope solitons only at lower 
or upper band edge, correspondingly. Thus long-lived 
BOs can be expected in the domain of parameters be- 
tween the two curves G^^'' (V) in FiglTJ 

Below we concentrate on the choice of the parame- 
ters corresponding to the specific points A, B, and C, 
in Fig. [T^. These are situations deviating form the opti- 
mally designed lattices, as is clearly seen in panels (b) and 
(c) of FiglH Indeed, while the points qM and are rel- 
atively close to each other, they do not exactly coincide, 
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FIG. 1: (a) Gj*' vs V which separate domains where 
x'^H^) ^ 0- The points A, B, and C correspond to the 
parameters explored below in Fig[2jD (point A), FigH^ and 
(point B), and FigOa (point C). Long living BO can be 
observed in the shadowed domain, (b) The inverse eflFec- 
tive mass AI~^ (solid line) and the effective nonlinearity x 
(dashed-dotted line) vs q. (c) Mx vs q for the parameters 
corresponding to the point B (panel a), qx ~ 0.5145 and 
qM ~ 0.5475 indicate the values of q, where x f^nd are 
respectively equal to zero. 



what results in the singularity of a{q) at the point qM 
and in the existence of the interval < q < qM where 
envelope solitons do not exist (however outside this in- 
terval (T is a relatively slow function of the wavevector). 

In Fig[2^ we present the time evolution of a small am- 
plitude envelope soliton obtained by direct numerical in- 
tegration of Eq. IH) . As initial condition we used the sta- 
tionary envelope soliton of Eq. ([1]) with 7 = 0, whose 
chemical potential belongs to the semi-infinite gap very 
close to the lowest allowed band. From this figure the 
existence of long-living BO with the period w 2 • 10'^ 
and spatial amplitude ~ 32.5tt, perfectly matching 
the semiclassical estimates 2/\^\ and (e'^'^'^ — e^~^)/(2|7|), 
respectively, is quite evident. Notice that the turning 
points of the spatial dynamics correspond to values of 
chemical potentials inside the gaps (very close to band 
edges) for which stationary solitons exist. In Fig[2t we 
have compared the shape of the soliton at the fifth right 
turning points of the BO with the profile of the sta- 
tionary state at the top of the band, from which we 
see that they practically coincide (the small difference 
is ascribed to the radiative effects of the BO dynam- 
ics). Remarkably accurate prediction of the theory is 
verified by comparing the dynamical profiles of the soli- 
ton at the turning points of the BO. The profiles at the 
times t = 10^; 3 • 10^; 5 • 10^; 7 • 10^ were practically in- 
distinguishable from the one depicted for t = 9 • lO'' in 
Figl2t. Numerical simulations performed on longer time 
scales (up to t = 2 • 10"*) showed no appreciable decay of 
the BO and perfect recovering of the soliton shape at the 
turning points. 

To appreciate the importance of the parameter design 




FIG. 2: (Color online) (a) Long-living BO of the envelope soli- 
ton for Gi — 1.0 (corresponding to point B in Fig[l^). (b) De- 
caying BO for Gi = 0.95 (corresponding to point A in Fig[l^). 
In both cases oscillations start near the bottom of the band 
at £ = —0.938 with parameter values 7 — —0.001, V = 3.0 
and with initial condition Xo = —65n. (c) The shape of the 
envelope soliton in panel (a) at time t = 9000 (thin solid lines) 
is compared with the stationary state (thick solid lines) with 
the same number of particles N = 2.42 and chemical potential 
nearby the upper edge of the band at e = —0.7324. (d) The 
chemical potential e{t) for long-living (Gi — 1.0, solid line) 
and decaying (Gi — 0.95, dashed line) BO. The lowest allowed 
band corresponds to the interval £ € [—0.93683, —0.73326]. 



implied by our theory, we depict in Figl^b the dynam- 
ics of a soliton corresponding to point A in FiglTJ We 
see that in this case the BO undergo fast decay due 
do spreading of the wave packet. The difference in the 
two types of BO is clearly seen also from the depen- 
dence of the chemical potential on time, computed as 
= iJZ.{\^.\'+U{x)\i;\'+g{x)\i;\')dx, with 

N = IV'l^'^a; denoting the number of atoms. This 
is shown in Fig[2]i from which we see that stable BO cor- 
respond to perfectly periodic trajectory (solid line) while 
decaying oscillations correspond to decays of the oscil- 
lations of the chemical potential (dashed line). Similar 
phenomena exist also if the BO is started from stationary 
states at the top, rather than at the bottom, of the band 
for the same parameter values of points B, C, in Fig[TJj,), 
as one can see from FiglSl 

In closing this letter we wish to discuss possible exper- 
imental settings for observing long-living BO in BECs 
trapped in uniformly accelerated linear and nonlinear 
OLs (with acceleration 27): L({X-jT'^) a,iidg{X~jT'^), 
respectively. The corresponding model equation in di- 
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FIG. 3: (Color online) Dynamics of BO for the cases (a) 
Gi = 1.0 and (b) Gi = 1.05 (corresponding, respectively, to 
points B and A in Fig[lji). In both cases the oscillations start 
near the top of the band at e = —0.7324, with 7 = —0.001, 
V = 3.0 and initial condition Xo = On. In (c) we compare 
the profile of the soliton at t = 9 ■ 10^ in panel (a) (thin solid 
lines) with the stationary solution (thick solid lines) with the 
same number of particles = 2.42 and belonging to the semi- 
infinite gap (e = —0.938). In (d) the dynamics of the chemical 
potential e{t) is shown for Gi = 1.0 (solid line) and Gi = 1.05 
(dashed line). 

mensionless variables takes the form 
i^fy = ~^xx+UiX - jT^)^' + g{x - 7r2)|*p^'. (4) 

After substituting V = e^'t'^-^^'^'^'^+'^'^'/^l^' and in- 
troducing new independent variables x — X — jT'^ and 
t = T, Eq. takes the form of a NLS equation with ex- 
ternal linear force ([1]). In this normalization the energy is 
measured in units of the recoil energy Er — fi^iT^ / {2m(P), 
where d is the lattice period and m is the mass of bosons, 
and the spatial and temporal variables are measured in 
units of (i/TT and h/Er, respectively. To check the exper- 
imental feasibility of the proposed setting, we consider 
a ^^Rb condensate in a trap with a transverse radial 
size a = 2 fim and with a period of linear and nonlin- 
ear lattices d ~ I /im. Then dimensionless parameters 
of long-living BO, depicted in Fig[2^ and [3^ will cor- 
respond to soliton, containing N « 3800 atoms, driven 
by external force 1.17 • lO^^^N (which is obtained, when 
the acceleration of linear and nonlinear lattices is of or- 
der of ~ 8mm/s). At the same time the above parame- 
ters of the nonlinear lattice could be created by the spa- 
tial variation (obtained by the optically induced Fesh- 
bach resonance) of the bosonic s-wave scattering length 
as{x) = ai°^ -I- ai^^ cos (i-nx/d) with the "average" value 



oh'^ = — 1.554 nm and the amplitude ai^^ = 2nm. For 
the parameters of FiglH the long-living BO will occur 
when the amplitude of the spatial variation of the scat- 
tering length is in the range 1.948 nm < a"P < 2.058 nm. 
This show that the long lived BO reported in this Let- 
ter can indeed be observed in the experimental settings 
available today. 
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